The Berry-Tabor conjecture for spin chains of Haldane-Shastry type 
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According to a long-standing conjecture of Berry and Tabor, the distribution of the spacings 
between consecutive levels of a "generic" integrable model should follow Poisson's law. In contrast, 
the spacings distribution of chaotic systems typically follows Wigner's law. An important exception 
to the Berry-Tabor conjecture is the integrable spin chain with long-range interactions introduced 
by Haldane and Shastry in 1988, whose spacings distribution is neither Poissonian nor of Wigner's 
type. In this letter we argue that the cumulative spacings distribution of this chain should follow 
the "square root of a logarithm" law recently proposed by us as a characteristic feature of all spin 
chains of Haldane-Shastry type. We also show in detail that the latter law is valid for the rational 
counterpart of the Haldane-Shastry chain introduced by Polychronakos. 

PACS numbers: 75.10.Pq, 05.45.Mt 



I. INTRODUCTION 

It is well known [1] that the distribution of (suitably 
normalized) spacings between consecutive eigenvalues for 
the Gaussian ensembles in random matrix theory is ap- 
proximately given by Wigner's surmise 

p{s) = (7r/2)sexp(-7rs^/4) . 

This behavior, which was first observed in the spectra 
of complex atomic nuclei, also seems to be characteristic 
of quantum chaotic systems like polygonal billiards 2]. 
On the other hand, for a "generic" quantum integrable 
system Berry and Tabor have conjectured [3] that the 
spacings distribution should instead follow Poisson's law 
p{s) = e"". Over the years, this conjecture has been 
verified for many integrable models of physical interest, 
such as the Heisenberg chain, the t-J model, the Hubbard 
model \M and the chiral Potts model Q. In a recent 
paper [6| it has been shown that there is an important 
exception to Berry and Tabor's conjecture, namely the 
integrable spin chain introduced by Haldane and Shastry 
in 1988 [3, lM]j whose spacings distribution follows neither 
Poisson's nor Wigner's law. It is the purpose of this 
letter to gain a deeper understanding of this fact, and to 
explore whether this property is shared by the analogous 
chain introduced by Polychronakos in ref . [9[ . 

The (antiferromagnetic) Haldane-Shastry (HS) spin 
chain describes a system of N spins equally spaced on 
a circle with long-range exchange interactions inversely 
proportional to the square of the chord distance between 
the spins. More precisely, its Hamiltonian is given by 
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where Sij is the operator permuting the i-th and j-th 
spins. (Unless otherwise stated, throughout the paper 
all sums and products run from 1 to A^.) The HS chain 
is closely related to the Hubbard model; for instance, the 
chain's ground state coincides exactly with Gutzwiller's 
variational wave function for the Hubbard model when 
the on-site interaction tends to infinity [l3, [HI, [l3 |. 
Another important characteristic of the HS chain is 
its connection with the Sutherland spin model of An 
type [13, Hj in, Ej [i3 ' from which it can be obtained 
by means of the so-called "freezing trick" [^ . This tech- 
nique essentially consists in taking the strong coupling 
limit in the spin Sutherland model, so that the internal 
and dynamical degrees of freedom decouple and the spins 
become frozen at the equilibrium positions of the scalar 
part of the potential, which are precisely the chain sites 
z9fc. This procedure can in fact be applied to all spin 
models of Calogero and Sutherland type, associated with 
both the Apf and BCn roots systems in Olshanetsky and 
Perelomov's schem e llSll. For instance, the spin Calogero 
model of Am type [13, [l3l 
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yields the Hamiltonian of the so-called Polychronakos- 
Frahm (PF) spin chain [I, HO] 



npF = Y.{^^-^,)-\l + s,,). 
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The sites of this chain are the coordinates of the unique 
critical point of the scalar part of the potential of the 
Hamiltonian ^ in the domain xi < • • • < xn, namely 
the numbers ^1 < • • • < ^at determined by the algebraic 
system 
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Remarkably, these numbers are just the roots of the Her- 
mite polynomial of degree N , as first pointed out by 
Calogero [2l|. 



In order to compute the spacings distributions of a 
spectrum, it is first necessary to transform the "raw" 
spectrum by applying the so-called unfolding map- 
ping [22l |. This mapping is defined by decomposing the 
cumulative level density F{£) as the sum of a fluctu- 
ating part F{\{£) and a continuous part 1]{£), which is 
then used to transform each energy £i, i = 1, . . . ,n, into 
an unfolded energy rji = ri{£i). In this way one obtains a 
uniformly distributed spectrum {r]i}f^i, regardless of the 
initial level density. One finally considers the normalized 
spacings Si = (77^+1 -?7i)/A, where A = {r]^ - m) / (n - 1) 
is the mean spacing of the unfolded energies, so that 
{si}^Zi has unit mean. The analysis of the spacings 
distribution has been carried out so far for the original 
(Ajv-type) HS chaiiij6|, its supersymmetric version j23l] 
and, very recently [24|, for the PF chain of BCN-tjpe 
introduced by Yamamoto and Tsuchiya 25]. The corre- 
sponding spacings distributions were found to be quali- 
tatively very similar, and differed essentially from both 
Wigner's and Poisson's distributions. 

For the PF chain of BCn type, we showed in ref. [2J| 
that for large N the cumulative spacings distribution 
P{s) = L p{t)dt was approximately given by 



P{s) c 1 



^_ /log(f^), (5) 



where Smax is the maximum spacing. In fact, the above 
approximation holds for any spectrum £1 = £niin < • • ■ < 
£n = i£'max obeying the following conditions: 
(i) The energies are equispaced, i.e., £i+i — £i — d for 

i = I, . . . ,n — 1. 
(ii) The level density (normalized to unity) is approxi- 
mately given by the Gaussian law 

5(^) = ^^e-^^, (6) 
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where /i and a are respectively the mean and the 
standard deviation of the spectrum. 

(iii) £max - Ai , M - ^min > CT- 

(iv) frnin and fmax a^c approximately symmetric with 
respect to ^, namely |£,nin + ■^max - 2/x| < £,„ax - 

As shown in ref. [2J], the above assumptions lead to 
the formula ^ with the following explicit expression for 
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It is relatively straightforward to check ^24i| that con- 
ditions (i)-(iv) are indeed satisfied by the PF chain of 
BCn type (with d = I). In the rest of this letter we shall 
discuss the applicability of the approximation ([5]) to the 
original (Ajy-typc) PF and HS chains. We shall see that 
conditions (i)-(iv) are all satisfied by the PF chain, so 
that eq. ([5]) is guaranteed to work in this case. On the 
other hand, we shall explain why the latter formula is 
still an excellent approximation for the HS chain, even if 
not all of the above conditions hold for this chain. 



II. THE POLYCHRONAKOS-FRAHM CHAIN 

The initial step in our analysis of the statistical prop- 
erties of the spectrum of the PF chain is the explicit 
knowledge of the partition function, which makes it pos- 
sible to compute the energy levels and their degeneracies 
for relatively large values of A^. As first shown by Poly- 
chronakos [26|], the partition function can be evaluated 
in closed form by applying the freezing trick to the spin 
Calogero model ^ . We shall present here an alternative 
expression for the partition function of the PF chain, 
which is more amenable to numerical computations than 
Polychronakos's original formula. 

The starting point in our derivation is the identity 



H ^W + 2anpF 
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where H'^'^ — H\ 

Calogero model. As explained in |2j|, this identity im- 
plies that for large a the energies of H are approximately 
of the form 
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Ef 
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where Ef^ and £j are two arbitrary eigenvalues of H^"^ and 
HpF- Although this relation cannot be used directly to 
compute the chain energies £j, one can use it to express 
the partition function Zpp of the PF chain as 



Zpf(T) = lim 



Z{2aT) 
Z^''{2aT) 
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in terms of the partition functions Z and Z'^'^ of H and 
H^'^, respectively. Both of these partition functions can 
be readily computed. Indeed, in the scalar case the en- 
ergies are given by p^ 



77'SC rp 



2a y rii 



Eo = aN{a{N- 1) + 1), (9) 



with n = (ni, . . . , Hat), where ni ^ • • • ^ un are nonneg- 
ative integers. Setting q = e~^/'^''^'^^ and neglecting the 
ground state energy Eq (which will also be subtracted 
from the energies of H) , we obtain 



Z"'=(2aT) = 



E' 
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, A^ — 1 , and pn = n^ , 



Y.Y.PJ ^ J2JPj^ we finally have 

E fi"-m^-^r'- (10) 

9i,...,Pjv^0 i 



The energies of the spin Hamiltonian H are still given by 
the RHS of ^, but now there is a degeneracy factor dn 
due to the spin [27] . More precisely, if m is the number 
of internal degrees of freedom {su{m) spin) and 
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Therefore 
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where Vn denotes the set of ordered partitions of N. 
Cahing again pj = Vj — Vj+i > 0, j = 1, . . . , r — 1, and 
Pr — Vr ^ 0, we have 

r r r r 

1=1 z— 1 j— i j — 1 

i 
where Xj = "^ki. Substituting the previous equation 

1=1 
into pip and proceeding as before we easily obtain 
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Since 1 ^ Ki < ■ ■ ■ < Kr-i < K^ = N, we can define 
N — r natural numbers K'^ by 

{K[,...,K'j^_^.}^{l,...,N-l}-{Ki,...,Kr-i}. 

From eqs. ^, (fTO|) and (fT2|) . it immediately follows that 
the partition function of the PF chain can be written as 

'' / \ 'V'^ K ^^^ 

2MT)=j:iiQ-^-' '•11(1-'^"') -(13) 

One can show that the latter expression is equivalent to 
Polychronakos's by arguing as in ref. [2^. Moreover, in 
the latter reference it is shown that eq. ([T^ implies that 
the energies of the PF chain are given by 
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where the motif 5 = {Si, . . . , Sn-i) is a sequence of O's 
and I's with at most to — 1 consecutive O's. From the 
previous formula it follows that the spectrum of the PF 
chain is a set of consecutive integers, so that condition 
(i) in the previous section is satisfied with spacing d — I. 
Hence the number of distinct energies n is simply the 
difference £max ~ i^min + 1, where the minimum and max- 
imum energies can be computed with the help of eq. pi)) . 
Indeed, the maximum energy is obviously obtained from 
the motif (1, . . . , 1), so that 



N-l 



fmax = 51 * = J ^(^ - 1) ■ 
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On the other hand, the minimum energy corresponds to 
the motif 



^=(...,1,0,...,0,1,0,...,0), 



(16) 



in which Si = 1 for i = N — jm with j = 1, . . . , fc = 
[N/m\. Thus 



fmin = J2^N - jm) ^kN~-k{k+ 1) 

J=l 
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We shall next examine the validity of the second con- 
dition in the previous section for the PF chain. In fact, 
eq. (I13p for the partition function turns out to be very ef- 
ficient for computing the spectrum of this chain for large 
values of N (up to A^ = 26 for to, = 2 or A^ = 20 for 
TO = 3, using Mathematica''"^ on a personal computer). 
In this way we have ascertained that the level density 
obeys the Gaussian law ([6]); cf. fig.[l]for the case iV = 25 
and TO = 2. The mean energy /i and the variance a in 




FIG. 1: Cumulative energy density (at its discontinuity 
points) and its continuous approximation J_ g{t)dt (con- 
tinuous red line) for the PF chain with TV = 25 and m — 2. 

eq. ([6]) can be computed in closed form for arbitrary val- 
ues of N and to, using several classical identities for the 
zeros of Hermite polynomials. More precisely, from the 



equality tr Si, 



TO 



N-l 



we easily obtain 
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where in the last equality we have used eq. (3.3b) of 
ref. [29|. In order to evaluate the last sum, we use the 
well-known identity ^29^, eq. (3.3a)] 
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Since, by antisymmetry, 



E 



e. - e, 



N{N - 1) 



multiplying eq. (J19p by ^^ and summing over i we obtain 



J2^^^1n{n^i)., 



and therefore 



4 V 771/ 

Similarly, using the identity 

and proceeding as in ref. [30|, we obtain 
2 _ tr(7^2^) (trHpp)' 



(20) 
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From [H, eq. (3.3d)] and ^ it follows that 



45$](C^-e.)-' = A^(iV-l)(2iV + 7)+^ef. (23) 
Multiplying eq. (J19p by ^^^ and summing over i we obtain 
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where we have used eq. (j20p and the obvious identity 
^^ ^, = 0. From eqs. @, (US)) and (glD we finally have 



.^ = ^(l-^)^(iV-l)(A^+^). (25) 

Equations dTS]), dTT]), (EH) and ^ imply that both 
(^max — /^)/o' and (/i — £niin)/cr gTow as A^^'^ whcu 
A^ — > oo, so that condition (iii) in the previous sec- 
tion is satisfied. The last of these conditions is also 
satisfied for large N, since by the equations just quoted 

l^min + fmax " 2^i\ = 0{N) whilc £„,ax " ^min = ^(iV^). 

By the discussion in the previous section, it follows that 
the cumulative spacings distribution of the PF chain is 
approximately given by eq. ([5]) for large N. We have in- 
deed verified that (O holds with great accuracy by com- 
puting P{s) for a wide range of values of N and m using 
eq. (I13p for the partition function. For instance, in the 
case N = 25 and to = 2 presented in fig. [2] the RHS 
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FIG. 2: Cumulative spacings distribution P{s) and its ana- 
lytic approximation ((5} (continuous red line) for the PF chain 
with N = 25 and m = 2. For convenience, we have also plot- 
ted Poisson's (green, long dashes) and Wigner's (green, short 
dashes) cumulative distributions. 



of (O fits the numerical data with a mean square error 
of 8.4 X 10-^ 

It should be noted that the parameter s,„ax in the for- 
mula ([5]) for P{s) can be computed explicitly as a func- 
tion of A^ and m using ([7]) (with d = 1), (PS)) . the identity 
n—l~ f max — f mill, and eqs. (fT5)) and p7|) . In particular, 
for large N the maximum spacing Smax behaves as 



V TO + 1 ^ ' 



just as for the PF chain of BCn type [2J. We thus 
conclude that for A'' ^ 1 the spacings distributions of 
the PF chains of An and BCn types are asymptotically 
equal, in spite of the fact that the spectra of these models 
are quite different. 



III. THE HALDANE-SHASTRY CHAIN 

It is well-known that the spectrum of the original HS 
chain ([T]) is not equispaced, so that the first condition 
in section 1 does not hold in this case. However, it was 
recently noted [2J| that the cumulative spacings distri- 
bution of this chain can be still approximated with great 
accuracy by a function of the form ([5]). The parameter 
Smax is given again by eq. ([7]), but now d is the spacing 
£i+i — £i with highest frequency and n is the total num- 
ber of energy levels. As an example, in fig. [3] we present 
the plot of P{s) and its approximation (O in the case 
A^ = 25 and to = 2, for which d — 2, n — 562, and hence 

Smax — O.llo. 

In the rest of this section we shall analyze why eq. ([5]) 
provides an excellent approximation to the cumulative 
spacings distribution of the HS chain, in spite of the fact 
that condition (i) in section 1 is not satisfied. We shall 
start by verifying that the remaining conditions (ii)-(iv) 
are satisfied in this case. In the first place, the fact that 




FIG. 3: Cumulative spacings distribution P{s) and its ana- 
lytic approximation © (continuous red line) for the HS chain 
with TV = 25 and m = 2. The mean square error of the fit is 
1.8 X 10"^ 



the level density is approximately Gaussian was estab- 
lished in ref. [6| , where it was also shown that the mean 
energy and its standard deviation are respectively given 

by 



12 V m/ 



a^ = —(l-^) N{N^ - 1)(N^ + 11) 
360 V m^J ^ A -r y 



(26) 
(27) 



As to the third condition, recall that the maximum en- 
ergy is given by [0| 



fmax = ^(A^'-l). 

6 



(28) 



The minimum energy can be computed from the analogue 
of eq. p^ , which for the HS chain reads [y, l3l| 



N-l 



£{S)^ J2S^^{N-^), 



i=l 

where again 5 is a motif with at most m ~ 1 consecutive 
O's. It can be shown that the motif yielding the minimum 
energy is again given by (fT6|) (or, alternatively, by the 
complementary motif S with (5^ — (5jv_i), so that 

k 
f min = ^ jm{N - jm) , 

with k = lN/m\ . Evaluating the latter sum we easily 
obtain 



£n 



6m 



{N ~1){N + m~l)iN -m + 2l), (29) 



where I = TV mod m. Equations ((26)) - ([29l) imply that (as 
for the PF chain) (Smax ^ m)/^ a-nd (/j, — i£niin)/<7 grow 
as iV^/^ when iV ^ oo, so that condition (iii) is clearly 
satisfied. Finally, from the previous equations it follows 



+ '?max - 2/x| = 0{N) while £,„ax - '?min = 

0{N^), which implies condition (iv). 

We shall next discuss in more detail to what extent 
the first condition fails. To this end, wc have used the 
formula for the partition function of the HS chain given 
in [y, eq. (22)] to compute the spectrum for a wide range 
of values of m and N such that rn^ ^ 2^^. Our results 
evidence that when N is sufficiently large the vast ma- 
jority of the differences di = £i+i — £i are equal to d = 1 
(for even N) ov d = 2 (for odd N), as shown in fig. ID 
Thus, we can say that the spectrum of the HS chain is 
quasi-equispaced for large N. Moreover, our computa- 
tions indicate that the differences di ^ d correspond to 
energies £i in the tail of the Gaussian distribution g{S); 
see, e.g., fig. [5] for the case m = 2, N = 26. 
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FIG. 4: Ratio ^^ ' , -* for the HS chain with m = 2, 3, 4, 5 



and N up to 28. 
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FIG. 5: Difi'erences di = £i+i — £i versus Si = {£i — iJ.)/cr for 
the HS chain with A'^ = 26 and m — 2. 

The two features of the spectrum of the HS chain just 
described provide the key to understanding why the for- 
mula ([5]) , with Smax given by ([7|) , still works remarkably 
well in this case. Indeed, since the level density is ap- 
proximately Gaussian, we have 



v{£) 



g{t)dt, 




FIG. 6: Comparison of the normalized spacings Si (blue dots) 
with the Gaussian distribution (|3Up (continuous red curve) for 
the HS chain with N = 26 and m — 2. 



Since the few values of di different from d correspond 
to energies £i several standard deviations away from /i, 
their associated spacings Si turn out to be very small. 
Therefore, we can write with great accuracy 



(n - l)d _(£i^^ 
s, ~ ^ , ^ e 2^'^ 



(gj-^-)^ 
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except for a few spacings Si very small compared to Smax; 
see, e.g., fig. |6]for the case A'' = 26, m = 2. As shown 
in ref. [2J|, eq. ([50)1 leads directly to the approximate 
formula ([5]) for the cumulative spacings distribution P{s). 
This explains why the latter formula is also valid for the 
HS chain, with the value of Smax given in eq. ([7]). 



so that 

By condition (iii) in section 1 we have rjn — 1 and 771 ~ 0, 
and thus A ~ l/(n— 1). Hence, the normalized spacings 
Si are approximately given by 



Si = 



m+i - Vi 



{n- l)dig{£i) 



{n~l)di _(£i-^o! 



2tt(t 



Acknowledgments 

This work was partially supported by the DGI under 
grant no. FIS2005-00752, and by the Complutense Uni- 
versity and the DGUI under grant no. GR74/07-910556. 
J.C.B. acknowledges the financial support of the Span- 
ish Ministry of Education and Science through an FPU 
scholarship. 



[1] T. Guhr, A. Miiller-Groeling, and H. A. Weidenmiiller, 

Phys. Rep. 299, 189 (1998). [18] 

[2] M. L. Mehta, Random Matrices (Elsevier, San Diego, 

2004), 3rd ed. [19] 

[3] M. V. Berry and M. Tabor, Proc. R. Soc. Lond. A 356, [20] 

375 (1977). [21] 

[4] D. Poilblanc, T. Ziman, J. Bellissard, F. Mila, and [22] 

J. Montambaux, Europhys. Lett. 22, 537 (1993). 
[5] J.-C. A. d'Auriac, J.-M. Maillard, and C. M. Viallet, J. [23] 

Phys. A 35, 4801 (2002). 
[6] F. Finkel and A. Gonzalez-Lopez, Phys. Rev. B 72, [24] 

174411(6) (2005). 
[7] F. D. M. Haldane, Phys. Rev. Lett. 60, 635 (1988). 
[8] B. S. Shastry, Phys. Rev. Lett. 60, 639 (1988). [25] 

[9] A. P. Polychronakos, Phys. Rev. Lett. 70, 2329 (1993). 
[10] J. Hubbard, Proc. Roy Soc. London Ser. A 276, 238 [26] 
(1963). [27] 

[11] M. C. Gutzwiller, Phys. Rev. Lett. 10, 159 (1963). 
[12] F. Gebhard and D. VoUhardt, Phys. Rev. Lett. 59, 1472 [28] 

(1987). 
[13] B. Sutherland, Phys. Rev. A 4, 2019 (1971). [29] 

[14] B. Sutherland, Phys. Rev. A 5, 1372 (1972). 
[15] Z. N. C. Ha and F. D. M. Haldane, Phys. Rev. B 46, 9359 [30] 

(1992). 
[16] K. Hikami and M. Wadati, J. Phys. Soc. Jpn. 62, 469 [31] 

(1993). 
[17] J. A. Minahan and A. P. Polychronakos, Phys. Lett. B 



302, 265 (1993). 

M. A. Olshanetsky and A. M. Perelomov, Phys. Rep. 94, 

313 (1983). 

F. Calogero, J. Math. Phys. 12, 419 (1971). 

H. Frahm, J. Phys. A 26, L473 (1993). 

F. Calogero, Lett. Nuovo Cimento 20, 251 (1977). 

F. Haake, Quantum Signatures of Chaos (Springer- 

Verlag, Berlin, 2001), 2nd ed. 

B. Basu-Mallick and N. Bondyopadhaya, Nucl. Phys. 

B757, 280 (2006). 

J. C. Barba, F. Finkel, A. Gonzalez-Lopez, and M. A. 

Rodriguez (2008), arXiv:0803.0922vl [cond-mat.stat- 

mech] . 

T. Yamamoto and O. Tsuchiya, J. Phys. A 29, 3977 

(1996). 

A. P. Polychronakos, Nucl. Phys. B419, 553 (1994). 

B. Basu-Malhck, H. Ujino, and M. Wadati, J. Phys. Soc. 
Jpn. 68, 3219 (1999). 

B. Basu-Mallick, N. Bondyopadhaya, K. Hikami, and 

D. Sen, Nucl. Phys. B782, 276 (2007). 

S. Ahmed, M. Bruschi, F. Calogero, M. A. Olshanetsky, 

and A. M. Perelomov, Nuovo Cimento B 49, 173 (1979). 

A. Enciso, F. Finkel, A. Gonzalez-Lopez, and M. A. 

Rodriguez, Nucl. Phys. B707, 553 (2005). 

F. D. M. Haldane, Z. N. C. Ha, J. C. Talstra, D. Bernard, 

and V. Pasquier, Phys. Rev. Lett. 69, 2021 (1992). 



